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In  this  paper,  transport  and  kinetic  parameters  of  lithium-ion  batteries  are  estimated  using  a  rigorous 
porous  electrode  theory  based  model.  The  rigorous  model  used  in  this  investigation  is  reformulated  using 
advanced  mathematical  techniques.  Since  batteries  and  other  electrochemical  devices  are  used  in  hybrid 
environments,  which  include  devices  with  time  constants  less  than  a  second  (like  supercapacitor),  we  need 
to  develop  parameter  estimation  codes  with  computation  time  less  than  a  second  or  a  few  milliseconds. 
In  this  investigation,  the  computation  time  for  parameter  estimation  measures  between  100  and  300  ms 
since  a  reformulated  battery  model  is  devised  especially  for  these  purposes.  Obtaining  the  numerical 
solution  for  battery  model  equations  is  very  difficult  towards  the  end  of  discharge  and  is  usually  neglected 
for  parameter  estimation  purposes.  However,  in  this  paper  the  estimation  takes  into  account  the  entire 
discharge  data  ranging  from  an  initial  potential  of  4.2  V  to  a  cut-off  potential  of  2.5  V.  It  is  found  from 
this  investigation  that  the  reformulated  lithium-ion  battery  model  is  efficient  and  accurate  in  estimating 
parameters. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Several  investigators  have  analyzed  various  rigorous  models  for 
lithium-ion  batteries  [1-12]  based  on  the  porous  electrode  theory 
coupled  with  concentration  solution  theory  and  modified  Ohm’s 
law.  There  are  some  important  advances  in  these  models  in  recent 
years  like  the  thermal  modeling  of  lithium-ion  batteries  [5]  and 
modeling  capacity  fade  mechanism  in  lithium-ion  batteries  [8].  The 
governing  equations  in  these  models  include  ten  nonlinear,  coupled 
and  multidimensional  partial  differential  equations  (in  the  spatial 
directions  x,  r  and  t)  that  are  needed  to  be  solved  simultaneously 
in  time  along  with  some  highly  nonlinear  algebraic  expressions 
for  transport  and  kinetic  parameters.  Rigorous  lithium-ion  battery 
models  need  any  where  from  several  seconds  to  minutes  to  sim¬ 
ulate  a  discharge  curve  depending  on  the  computer,  solver,  etc. 
This  computational  difficulty  in  using  rigorous  battery  models  is 
due  to  the  large  number  of  equations  that  result  from  finite  dif¬ 
ference  reformulation  of  battery  models.  Table  1  shows  the  total 
number  of  differential-algebraic  equations  (DAEs)  needed  for  the 
rigorous  model  and  the  computation  time  required  when  the  model 
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equations  are  solved  for  one  charge  or  discharge  curve  and  for  one 
complete  cycle.  All  the  simulations  in  this  work  are  based  on  a 
PC  with  2.5  GHz  processor  and  2  GB  RAM.  Rigorous  version  of  the 
model  is  not  ideal  due  to  the  requirement  of  sufficient  computa¬ 
tional  resources.  Whereas,  the  reformulated  efficient  model  could 
eventually  in  future  run  in  a  microprocessor  to  estimate  parameters 
in  a  few  milliseconds. 

Until  recently  [13],  there  were  no  significant  efforts  for  esti¬ 
mating  parameters  using  the  rigorous  battery  model  because 
of  computational  constraints.  There  are  several  mathematical 
investigations  for  rapid  solution  of  complicated  model  equations 
[14-22].  A  more  detailed  article  that  presents  the  complete  mathe¬ 
matical  theory  for  an  effective  approach  to  simulate  physics-based 
lithium-ion  battery  models  in  real-time  (milliseconds)  is  under 
the  process  of  publication  [23].  Parameter  estimation  helps  in 
calculating  capacity  fade  and  life  time  of  the  battery.  However, 
empirical  models  or  single  particle  model  cannot  help  in  estimating 
accurate  parameter  values.  The  capacity  fade  or  the  state-of-health 
information  obtained  using  these  parameters  will  not  guarantee 
successful  prediction  for  applications  like  hybrid  vehicles,  satellites 
or  military  applications.  Also,  it  is  to  be  noted  that  the  batteries 
and  other  electrochemical  devices  are  now-a-days  emerging  as  a 
promising  candidate  for  hybrid  power  systems.  Since  in  a  hybrid 
environment,  there  are  devices  like  the  supercapacitor  with  time 
constant  less  than  a  second,  we  also  need  to  have  a  battery  model 
that  can  be  simulated  in  few  milliseconds.  To  this  point,  this  inves¬ 
tigation  and  our  future  efforts  on  lithium-ion  battery  models  are 
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Nomenclature 

specific  surface  area  of  electrode  i  (i=p,  n)  (m2  m-3) 

brugg, 

Bruggman  coefficient  of  region  i  ( i  =  p ,  s,  n) 

c 

electrolyte  concentration  (molm-3) 

Co 

initial  electrolyte  concentration  (molm-3) 

Cs  ,i 

concentration  of  lithium  ions  in  the  intercalation 
particle  of  electrode  i  (i=p,  n)  (molm-3) 

Cs.i.O 

initial  concentration  of  lithium  ions  in  the  interca¬ 
lation  particle  of  electrode  i  (z=p,  n)  (molm-3) 

Cs.i.max 

maximum  concentration  of  lithium  ions  in  the  inter¬ 
calation  particle  of  electrode  i  (z=p,  n)  (molm-3) 

D 

electrolyte  diffusion  coefficient  (m2  s-1 ) 

Ds,i 

lithium  ion  diffusion  coefficient  in  the  intercalation 
particle  of  electrode  i  (i =p,  ri)  (m2  s-1 ) 

F 

Faraday’s  constant  (C  mol-1 ) 

h 

solid-phase  current  density  (Am-2) 

h 

solution  phase  current  density  (Am-2) 

h,0 

exchange  current  density  for  the  solvent  reduction 
reaction  (Am-2) 

I 

applied  current  density  (A cm-2) 

it 

wall  flux  of  Li+  on  the  intercalation  particle  of  elec¬ 
trode  i  ( i  =  n,  p)  (mol  m-2  s-1 ) 

is 

solvent  reduction  current  density  (mol  m-2  s-1 ) 

intercalation/deintercalation  reaction  rate  constant 
of  electrode  i  (i=p,  n)  (mol  (molm-3)-15) 

/, 

thickness  of  region  i  (i =p,  s,  n)  (m) 

Ms 

molar  weight  of  the  solvent  reaction  product 
(gmol-1) 

n 

negative  electrode 

P 

positive  electrode 

r 

radial  coordinate  (m) 

R 

universal  gas  constant  (J  (mol-1  K-1 )) 

Ri 

radius  of  the  intercalation  particle  of  electrode  i 
(i=p,n)(  m) 

Rsei 

initial  SEI  layer  resistance  at  the  negative  electrode 

(£2  m-2) 

s 

separator 

t+ 

Li+  transference  number  in  the  electrolyte 

T 

absolute  temperature  (I<) 

Ui 

open  circuit  potential  of  electrode  i  (i =p,  n)  (V) 

Us 

standard  potential  of  the  solvent  reduction  reaction 

(V) 

X 

spatial  coordinate  (m) 

Greek  symbols 

<pi 

solid-phase  potential  (V) 

@2 

electrolyte  phase  potential  (V) 

8 

thickness  of  the  solvent  reduction  product  film  (m) 

initial  thickness  of  the  solvent  reduction  product 
film  (m) 

£i 

porosity  of  region  i  (i =p,  s,  n) 

Sf ,i 

volume  fraction  of  fillers  of  electrode  i  (i=p,  ri) 

K 

ionic  conductivity  of  the  electrolyte  (S  m-1 ) 

effective  ionic  conductivity  of  the  electrolyte  in 
region  i  (i =p,  s,  n)  (S  m-1 ) 

dimensionless  concentration  of  lithium  ions  in  the 
intercalation  particle  of  electrode  i  (0,-  =  cStj/cSf,-fmax) 

Ps 

density  of  the  solvent  reduction  product  film 

(gm~3) 

electronic  conductivity  of  the  solid  phase  of  elec¬ 
trode  i  (i=p,  n)  (Sm-1) 

effective  electronic  conductivity  of  the  solid  phase 
of  electrode  i  (i=p,  n)  (Sm-1) 

Table  1 

Properties  of  the  rigorous  porous  electrode  model  and  the  reformulated  porous 
electrode  model 


Computation  details 

Rigorous  model 

Reformulated  model 

Number  of  equations  for 
charge  or  discharge  curve 

4800  DAEs 

49  DAEs 

Computation  time  for  charge 
or  discharge  curve 

1-2  min 

15-45  ms 

Estimated  time  for  one 
complete  cycle 

5-10  min 

75-100  ms 

focused  on  real-time  simulation,  parameter  estimation,  dynamic 
optimization,  and  on-line  control  and  monitoring  of  batteries  that 
can  work  with  CPU  time  less  than  a  second. 

Santhanagopalan  et  al.  [13]  used  the  Levenberg-Marquardt 
method  to  estimate  parameters  using  rigorous  lithium-ion  battery 
models.  Two  different  models,  the  single  particle  model  and  the  rig¬ 
orous  porous  electrode  model,  were  used  to  estimate  parameters  at 
four  different  rates.  Confidence  intervals  for  the  estimated  param¬ 
eters  were  reported,  however,  the  computational  time  involved 
in  this  investigation  was  not  highlighted;  perhaps  because  of  the 
necessity  to  do  expensive,  cumbersome  and  time-consuming  simu¬ 
lation  while  using  the  porous  electrode  model.  The  authors  [19,20] 
have  studied  the  validity  of  estimating  the  solid-phase  diffusion 
coefficient  of  a  lithium  intercalation  electrode  from  impedance 
measurement  by  a  modified  electrochemical  impedance  spec¬ 
troscopy  method.  This  paper  showed  that  single  particle  models 
are  not  accurate  and  are  hence  not  useful. 

To  overcome  this  drawback,  we  have  devised  a  novel  reformu¬ 
lated  lithium-ion  battery  model  that  can  help  in  simulating  battery 
behavior  in  real-time  (milliseconds)  [21].  This  novel  reformulated 
battery  model  has  been  used  in  this  article  for  parameter  estima¬ 
tion.  The  model  used  in  this  investigation  is  based  on  the  same 
porous  electrode  model  for  the  lithium-ion  cell  that  is  used  in  ref. 
[13].  However,  it  is  found  that  the  computation  time  required  for 
the  parameter  estimation  is  only  on  the  order  of  a  few  millisec¬ 
onds  in  this  investigation.  Table  1  also  shows  the  number  of  DAEs 
required  in  solving  the  reformulated  battery  model  and  the  com¬ 
putation  time  associated  with  it  while  solving  the  model  for  one 
charge  or  discharge  curve  and  for  one  complete  cycle.  The  refor¬ 
mulated  lithium-ion  battery  model  is  efficient  and  accurate  for 
parameter  estimation.  It  is  also  important  to  note  that  with  the  use 
of  reformulated  lithium-ion  battery  model  the  numerical  Jacobian 
calculation  is  good  enough  for  accurate  parameter  estimation.  Thus, 
this  investigation  enables  the  possibility  for  real-time  estimation  of 
parameters  that  can  be  useful  in  predicting  the  capacity  fade  and  life 
time  of  the  lithium-ion  batteries.  It  should  be  noted  that  the  model 
reformulation  provides  predictive  models  and  are  as  good  as  origi¬ 
nal  model  for  a  wide  range  of  operating  conditions.  It  is  significantly 
different  from  proper  orthogonal  decomposition  approach  which 
is  based  on  curve  fitting  eigen  values  and  eigen  vectors  to  produce 
reduced  order  models  from  rigorous  numerical  simulation  [24-26]. 

2.  Reformulated  lithium-ion  battery  model 

The  existing  rigorous  physics-based  lithium-ion  battery  model 
[5]  has  multiple  partial  differential  equations  in  multiple  domains. 
The  lithium-ion  battery  model  involves  two  independent  spatial 
variables  (radial  direction,  r  and  linear  direction,  x)  and  a  time 
variable,  t.  There  are  four  dependent  variables,  namely  electrolyte 
concentration  (c),  electrolyte  potential  (#2).  solid-phase  concen¬ 
tration  (cs)  and  solid-phase  potential  (<2>i)  for  each  of  the  two 
porous  electrodes.  The  two  dependent  variables  for  separator  are 
electrolyte  concentration  (c)  and  electrolyte  potential  (<3>2)-  The 
governing  equations  for  all  these  variables  are  shown  in  Table  2 
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Table  2 

Governing  equations  for  a  lithium-ion  cell  based  on  the  porous  electrode  theory 
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Region  Eq.  no. 


SP  —  ^eff,p  §^§  +  Qp(l  -  t+l/p  i 

rr  rr  /r  rr  I  2Xeff+RG  i  /-  'j  3  Inc  _  r 

aeff,p  9x  *eff,p  3X  +  F  D  f+J  dx 

°eff,p  —  apfjp 

3cs 

at 


Cathode  (1) 
(2) 

(3) 

(4) 

Separator  (5) 
(6) 

Anode  (7) 
(8) 

(9) 

(10) 


p  initial  condition  i 

n  RT, 


eff,p  ~faT  —  uPrJP 

W  =  {r2w)  initial  condition  cs  1 1=0 

P  3C  _  n  rr  3^C 

Ss  at  _  Ueff+  dx2 


J  —  ee  302  I  2yeff,sRT/-|  t  >,3  Inc 

^eff,s  9x  +  F  D  r+J  dx 
Sn  §  =  Defffn  0  +  an(l  -  4)/n  initial  condition  ci 

90i  f/.  302  ,  2/ceffnRT  M  y-iainc  i 

_oreff,n  -5T  “  *eff,n  -3*-  +  “ 1 F - ( 1  _  )  ~dx~  = 


t=0  =  Co 


4ff,p  a9x2  lx=0  -  0  &  4ff,f 


301  I  _  J  o  _  30i  I  _  n 

dx  |x=0  "  ^eff.p  eff’p  3x  |x=/p  “ 

=  °-5cs,max,p  ^  |r=Q  =  0  &  JP=  Ds>P  jjf  |r=Rp 

— ^eff,p|^|  .  =— ^eff,s|^|  .  &  —  Oeff,s  ^  I  ,  , 

’  ^  I  x=/p ox  lx=/p,+  ax  lx=/p+/s,- 

302  I  o,  ^  302  I 

1  ®  -  /ceff  s 

’  dx  lx=/p+/s,- 


-  Deff,n  9x 


Keff-p  8S2  L=ip,- _  &2  L=iP,+ 

=-D^"lL=,P+,s,+&-D^i|teIp+ls+l„ 

=  0 


X=/p+/s ,+ 

m-K*  ^1 

-  «efr,n  3x  |x=(p+Js  . 


^eff,!!  9x21  —  ®nTjn 
+  =  4  I  (r2  +  )  initial  condition  cs|t=o  =  0.85cs,max,n  f^|r=0 


=-^,n^L=(p+,s+&tP2|x=Jp+,s+,n 


=  0 


-/Ceff  s 

ett’S  ^x  |x=/p+/j 

30i  I  _  n  o  30i  I 

eff.n  dx  L=lp+Is  “  Hr.  |x=/p+/s+ln 

=  0&Jn  =  -DS,„f  |r=Rn 


_ / 

°eff,n 


Table  3 

Expressions  used  in  the  lithium-ion  battery  model  given  by  Table  1 


Keffti  =  ebruggi(4  1253  xl(T2+5.007  x  1CT4c  -  4.7212  x  1CT7c2  +  1.5094  x  l(T10c3  -  1.6018  x  10“14c4),  i  =  p,s,n 

<re ff,|  =  ciO  -  Ci  -  £ffi)t  i  =  P,n 

Deff  j  =  DF^ru§§> ,  i~p,s,n 

’  3  *  x  • 

Of  =  ^-(1  -Si-Efi),  i=p,n 


0.5  F 

-(0^-02  -Up) 


jp  =  2/cp(cs,max,p  -  Cs,p|r=i?p)0'5Cs,p|r=j?p0-5C0-5  sin/i 

-4.656  +  88.6696>2  -401. 1196>4  +  342.909$;  -462.471$  +  433.434$° 
U  _  _ v _ v _ v _ p _ p_ 

P  ~~  -1.0  +  18.93302  -  79.5326>4  +  37.3116>°  -  73.0836>3  +  95.966>J° 

where  :  0P  =  cSlP  I  /cS)P  ,max 
I  r=i?n 


j n  =  2/cn(Cs,max,n  _  Cs  n  )  Cs  n 
I  r=Rn  I 


c0-5  sin  h 


0.5  F 
RT 


(0\  -  02  ~  Un  +  ERsEt/n ) 


Un  =  0.7222  +  0.13876>n  +  0.029^-5  -  °'0'172  +  Q-0^19  +  0.2808  exp(0.90  -  156>n)  -  0.7984  exp(0.44656>n  -  0.4108) 

On  ^n1'5 


where  :  0n  =  cs,n  „  /cs,n,max 
I  r=i?n 


and  the  expressions  used  for  the  electrochemical  reaction  kinetics 
and  other  transport  properties  can  be  found  in  Table  3.  A  review 
on  available  simulation  techniques  for  solving  such  a  complicated 
lithium-ion  battery  model  was  made  earlier  [9]. 

These  model  equations  are  usually  discretized  using  finite  differ¬ 
ences  in  the  spatial  coordinates  resulting  in  4800  DAEs.  When  4800 
DAEs  are  needed  to  solve  for  a  simulation  of  the  model  equations, 
the  Jacobian  for  parameter  estimation  will  add  some  thousands 
more  DAEs.  Together,  we  will  have  several  thousands  of  DAEs  to 
estimate  one  parameter.  If  we  go  for  instead  to  predict  more  than 
one  parameter,  then  the  number  of  equations  to  be  solved  will 
keep  increasing  proportionally.  Obviously,  this  situation  is  not  ideal 
for  emerging  applications  like  hybrid  power  systems,  on-line  con¬ 
trol  and  monitoring  of  batteries  and  other  electrochemical  power 
sources.  It  will  take  several  hours  or  days  to  estimate  parameters  in 
the  model  equations. 

The  authors  have  worked  extensively  in  model  reformulation 
and  have  already  published  the  details  on  the  reformulation  of 
lithium-ion  battery  model  [21,22].  This  can  enable  a  simulation 
of  porous  model  equations  in  15-45  ms  as  against  90-120  s  for  a 
similar  simulation  using  rigorous  porous  model  equations.  This 
is  achieved  by  efficiently  combining  the  advanced  mathematical 
approximation  techniques  with  the  actual  physics  of  the  electro¬ 
chemical  system  under  consideration.  Now,  the  number  of  DAEs 
needed  for  a  simulation  is  only  29-49. 


3.  Parameter  estimation 

A  list  of  fixed  parameters  for  the  model  equations  are  given 
in  Table  4,  and  the  following  parameters  are  estimated  (1)  the 
solid-phase  diffusion  coefficient  in  negative  electrode,  Ds,n  (2)  the 
electrolyte  diffusion  coefficient,  D  (3)  the  electrochemical  rate  con¬ 
stant  in  positive  electrode,  kp  and  (4)  the  electrochemical  rate 
constant  in  negative  electrode,  kn.  The  unknown  parameters  are 
estimated  using  Gauss-Newton  method  from  the  synthetic  experi¬ 
mental  data  obtained  using  computer  generated  random  numbers. 
Here,  the  synthetic  data  is  obtained  by  distributing  5%  randomness 
error  to  the  actual  model  predicted  discharge  data  by  consider¬ 
ing  102  data  points.  The  expression  used  to  generate  experimental 
values  using  Maple’s  random  number  generator  is, 

-^experimental  =  ^analytical  [  0 .95  +  0.01  (Random  number 
between  0  and  10)]. 

The  synthetic  experimental  data  are  then  provided  to  the  param¬ 
eter  estimation  algorithm  along  with  appropriate  initial  guesses. 
The  following  steps  are  involved  in  the  estimation  of  parameters: 
(i)  start  with  good  initial  guesses  for  parameters,  k(0)  (ii)  compute 
the  state  and  sensitivity  variables  of  the  model  at  each  data  point 
and  set  up  the  vector  with  experimental  values  Yexp,  predicted 
values  Ypre  and  the  Jacobian  matrix  J  (iii)  the  correction  fac- 


364 


V.  Boovaragavan  et  al.  /  Journal  of  Power  Sources  183  (2008)  361-365 


Table  4 

Parameters  used  for  the  simulation  (LiCo02  and  LiCg  system) 


Symbol 

Unit 

Positive  electrode 

Separator 

Negative  electrode 

Sirr1 

100 

100 

ef ,i 

0.025 

0.0326 

0.385 

0.724 

0.485 

Brugg 

4 

Ds  ,i 

m2  s-1 

l.OxlO-14 

3.9  x  10“14 

D 

m2  s-1 

7.5  x  10-10 

ki 

mol(s-1  m-2)[(molm_3)1+aa'1]-1 

2.334  x  10-11 

5.0307  x  10-11 

cs,i,max 

mol  ITT3 

51,554 

30,555 

cs,i0 

mol  ITT3 

0.4955  x  51,554 

0.8551  x  30,555 

Co 

mol  ITT3 

1000 

RP 

m 

2.0  x  10-6 

2.0  x  10-6 

k 

m 

00 

o 

X 

o 

G) 

25  x 10-6 

00 

00 

x 

o 

CT5 

Rsei 

Q  m2 

0.0 

T+ 

0.363 

F 

Cmol-1 

96,487 

R 

J  (mol-1  K-1) 
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tor  is  obtained  by  using  the  expression,  Ak  =  (JTJ)_1JT(YeXp  -  Ypre) 
(iv)  using  this,  the  updated  parameter  values  can  be  obtained  as 
k(l+1)  =  k(l)  +  rAk(l),  the  predicted  parameter  for  the  next  iteration 
is  r  is  the  correction  factor.  Steps  (ii)  to  (iv)  are  repeated 

until  a  required  accuracy  is  reached  or  the  difference  between  the 
parameter  value  of  current  iteration  and  previous  iteration  is  zero. 

In  this  paper,  the  95%  confidence  interval  is  calculated  for  the 
transport  and  kinetic  parameters  by  using  the  following  expression 
[27,28].  p±tN_n  a/2{s/VN  -  n)y/aji  where  p  is  the  point  estimate 
for  the  parameter  p,  t  is  a  value  of  Student’s  t  distribution  with 
(N  -  n)  degrees  of  freedom  and  a  =  0.5,  i.e.,  95%  confidence,  N  is  the 
total  number  of  discharge  data  points  and  n  is  the  number  of  fitting 
parameters,  an  is  the  ith  element  of  the  principal  diagonal  of  (JTJ)_1 , 
and  s  is  variance  or  deviation  variable  that  can  be  calculated  by  s2  = 
jfjLt  [('''cell);  -  (V'ceii);*]2  where  (VceiOj  and  (Vcen)*  are  the  predicted 
and  experimental  values  of  the  discharge  curve,  respectively,  at  the 
jth  point  at  time  t. 

4.  Results  and  discussion 

The  parameters  were  estimated  for  the  following  three  different 
cases:  (1 )  the  solid-phase  diffusion  coefficient  alone— estimation  of 
one  parameter;  (2)  the  solid-phase  diffusion  coefficient  and  elec¬ 
trolyte  diffusion  coefficient— estimation  of  two  parameters;  (3)  the 
solid-phase  diffusion  coefficient,  electrochemical  rate  constant  at 
the  positive  electrode  and  electrochemical  rate  constant  at  the  neg¬ 
ative  electrode— estimation  of  three  parameters. 

4.1.  One  parameter 

Parameter  estimation  was  carried  out  to  compare  the  compu¬ 
tation  statistics  for  the  rigorous  model  and  efficient  reformulated 
model.  The  parameter  estimated  was  the  solid-phase  diffusion 
coefficient  in  the  cathode;  its  true  value  is  1  x  10-15.  Table  5 
shows  that  the  rigorous  model  has  10  PDEs  or  4800  DAEs  and 
takes  1  -2  min  to  complete  a  single  iteration.  It  took  4  iterations 
to  converge  (using  a  very  good  initial  guess)  and  the  total  time  for 
estimating  one  parameter  is  5-10  min  using  an  efficient  solver/good 
initialization.  The  reformulated  model  has  49  DAE/IDEs  and  takes 
only  60  ms  to  estimate  a  parameter  in  FORTRAN  environment.  The 
manual  time  taken  between  running  the  codes  at  different  updated 
parameter  values  is  ignored  here  for  the  purpose  of  stable  esti¬ 
mation  algorithm.  The  synthetic  data  (for  1C  discharge)  used  are 
shown  in  Fig.  1  in  comparison  with  simulated  model  data  and  the 
initial  guess. 


4.2.  Two  parameters 

The  computational  run  time  was  determined  for  estimating  two 
parameters  in  the  model  equations.  Fig.  2  represents  the  discharge 
curves  for  the  synthetic  experimental  data,  the  discharge  curve 
based  on  the  initial  guesses  and  the  converged  discharge  data. 
Table  5  shows  the  computation  time  taken  for  the  estimation  of 
two  parameters. 

4.3.  Three  parameters 

The  computational  run  time  was  determined  for  estimating 
three  parameters  in  the  model  equations.  Fig.  3  represents  the  dis- 


Table  5 

Estimation  of  solid-phase  diffusion  coefficient,  electrolyte  diffusion  coefficient  and 
electrochemical  reaction  rate  constants  in  the  positive  and  negative  electrodes 


Number  of  parameters 

Runtime  for  each  iteration 

Number  of 
iterations 

One  parameter  estimation 

80s  (rigorous  model) 

15  ms  (reformulated  model) 

4 

Two  parameters  estimation 

32  ms  (reformulated  model) 

5 

Three  parameters  estimation 

35  ms  (reformulated  model) 

2 

Comparison  of  computational  time  taken  by  the  rigorous  porous  electrode  model 
and  the  reformulated  porous  electrode  model. 


Fig.  1.  Estimation  of  solid-phase  diffusion  coefficient  alone.  Comparison  of  synthetic 
experimental  data,  model  predicted  discharge  data  based  on  the  initial  guess  and 
model  converged  discharge  data. 
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Fig.  2.  Estimation  of  solid-phase  diffusion  coefficient  and  electrolyte  diffusion  coef¬ 
ficient.  Comparison  of  synthetic  experimental  data,  model  predicted  discharge  data 
based  on  the  initial  guess  and  model  converged  discharge  data. 


Fig.  3.  Estimation  of  solid-phase  diffusion  coefficient  and  electrochemical  rate 
constants  in  the  positive  and  the  negative  electrodes.  Comparison  of  synthetic  exper¬ 
imental  data,  model  predicted  discharge  data  based  on  the  initial  guess  and  model 
converged  discharge  data. 


Table  6 

Confidence  intervals  for  the  point  estimation  of  transport  and  kinetic  parameters  of 
the  lithium-ion  battery 


charge  curves  for  the  synthetic  experimental  data,  the  discharge 
curve  based  on  the  initial  guesses  and  the  converged  discharge  data. 
Table  5  shows  the  computation  time  taken  for  the  estimation  of 
three  parameters. 

The  above  three  cases  show  that  the  reformulated  lithium-ion 
battery  model  is  very  powerful  and  computationally  efficient  for 
parameter  estimation.  Table  6  shows  the  confidence  interval  for 
the  transport  and  kinetic  parameters  that  are  estimated  using  this 
efficient  reformulated  lithium-ion  battery  model. 


5.  Conclusions 

This  investigation  reports  a  real-time  parameter  estimation 
using  a  reformulated  lithium-ion  battery  model.  The  computational 
benefits  involved  in  this  approach  are  several  times  better  than 
the  current  state-of-the-art.  These  results  indicate  that  the  rigor¬ 
ous  optimization  for  optimal  operating  conditions  of  lithium-ion 
battery  model  can  now  be  evaluated  using  this  reformulated  model 
[29-31  ].  It  is  to  be  noted  in  this  work  that  the  estimation  of  param¬ 
eters  using  lithium-ion  battery  models  are  performed  only  for  up  to 
2C  rate  of  discharge.  Further  reformulation  of  the  rigorous  battery 
model  is  required  to  enable  estimation  of  parameters  at  high  rates 
of  discharge.  Work  in  this  regard  is  currently  in  progress.  Future 
work  will  also  involve  optimization  of  lithium-ion  battery  model 
and  other  aspects  related  to  the  control  of  electrochemical  devices. 
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